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Recent numerical studies indicate that uncertainties in the treatment of baryonic physics can affect 
predictions for weak lensing shear power spectra at a level that is significant for several forthcoming 
surveys such as the Dark Energy Survey (DES), the SuperNova/Acceleration Probe (SNAP), and 
the Large Synoptic Survey Telescope (LSST). Correspondingly, we show that baryonic effects can 
significantly bias dark energy parameter measurements. Elimination of such potential biases by 
neglecting information in multipoles beyond several hundred leads to weaker parameter constraints 
by a factor of ~ 2 — 3 compared with using information out to multipoles of several thousand. 
Fortunately, the same numerical studies that explore the influence of baryons indicate that they 
primarily affect power spectra by altering halo structure through the relation between halo mass 
and mean effective halo concentration. We explore the ability of future weak lensing surveys to 
constrain both the internal structures of halos and the properties of the dark energy simultaneously 
as a first step toward self calibrating for the physics of baryons. In this approach, parameter biases 
are greatly reduced and no parameter constraint is degraded by more than ~ 40% in the case of LSST 
or 30% in the cases of SNAP or DES. Modest prior knowledge of the halo concentration relation and 
its redshift evolution greatly improves even these forecasts. In addition, we find that these surveys 
can constrain effective halo concentrations themselves usefully with shear power spectra alone. In 
the most restrictive case of a power-law relation for halo concentration as a function of mass and 
redshift, the concentrations of halos of mass m ~ 10" /i-^Mq at z ~ 0.2 can be constrained to 
better than 10%. Our results suggest that inferring dark energy parameters through shear spectra 
can be made robust to baryonic physics and that this procedure may even provide useful constraints 
on galaxy formation models. 

PACS numbers: 98.80.-k,98.62.Py,98.35.Gi 



I. INTRODUCTION 

An ever-expanding set of observational data indicate 
that roughly ~ 75% of the energy in the universe is in 
the form of dark energy whose negative pressure drives 
an accelerated cosmological expansion (e.g., [U, 0, 0, 0, H, 
0)0,11, Hi)- Determining the nature of the dark energy is 
one of the most profound problems facing physicists and 
astronomers. Numerous contemporary and forthcoming 
experiments aim to shed light on the dark energy problem 
using a variety of techniques. One of the most promis- 
ing probes of dark energy is weak gravitational lens- 
ing. Forthcoming weak lensing surveys expect to measure 
the statistics of the matter density fluctuation field with 
exquisite precision. Dark energy affects both the rate of 
structure growth and the relative angular diameter dis- 
tances between different redshifts and so these precise 
determinations of the matter fluctuation spectrum bring 
the promise of strin gent constraints on dark energy pa- 
rameters [IS M MM M M, M M, M, M, Mm ■ 

Of course, exploiting precise measurements of fluctua- 
tions in the density field to constrain dark energy requires 



precise predictions for these inhomogeneities. Maximiz- 
ing the constraining power of future surveys in this re- 
gard demands that predictions for the linear and non- 
linear matter power spectra be accurate at the percent 
level [l^, [2^. This remains a challenge for numerical 
simulations even in the absence of baryonic physics [2^ , 
but this challenge should yield to unrelenting increases 
in computational power. Perhaps more important is the 
treatment of the baryonic component of the universe. 
Early analytic studies indicated that baryonic processes 
may have a significant influence on shear power spectra 
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25,l26|. More recently, Jing et al. [27| and Rudd et al. 

28| showed that numerical simulations that include bary- 
onic processes predict matter power spectra that differ 
from those that result from dissipationless iV-body cal- 
culations at a level that exceeds the precision of various 
forthcoming experiments. The evolution of the baryonic 
component of the universe cannot be modeled directly 
with current computational limitations and this circum- 
stance is unlikely to yield to increases in computational 
power in the foreseeable future. As a consequence, all 
calculations rely on effective models that attempt to ap- 
proximate the net, large-scale influences of numerous pro- 
cesses that occur on scales far below the numerical res- 
olution that may be attained with any simulation. The 
implication of Refs. [27], |28i] is that the ability of weak 
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lensing surveys to constrain dark energy parameters may 
be severely limited by this inability to account for bary- 
onic processes reliably. 

The goal of this study is to investigate this implica- 
tion. In particular, we aim to determine if measurements 
of shear power spectra (including tomography with pho- 
tometric redshifts) contain sufficient information to con- 
strain both the influence of baryons and the properties of 
the dark energy. Rudd et al. [28] found that the influence 
of baryons is due almost entirely to modiflcations in the 
matter distribution within dark matter halos. We model 
the influence of baryons in this way, so our calculation is 
tantamount to determining the ability of future surveys 
to constrain a physically-motivated model for the density 
profiles of dark matter halos and dark energy parameters 
simultaneously. Naturally, gravitational lensing is sensi- 
tive to all matter along the lines of sight to sources. In the 
following we refer to the distribution of mass within halos 
for simplicity, but we do not intend to restrict our con- 
sideration to the dark matter. We use this terminology 
as a shorthand to refer to the total mass distribution of 
the composite dark matter and baryonic systems within 
halo virial radii. 

Our results suggest that forthcoming surveys will have 
sufficient information to calibrate for the influence of 
baryons on the matter power spectrum using weak lens- 
ing shear power spectra alone, with minimal degrada- 
tion in the constraints on dark energy parameters. In 
fact, we study models that are quite general so our re- 
sults suggest that such surveys can also calibrate ad- 
ditional systematics that may afflict a similar range of 
scales. Of course, one need not regard halo structure 
strictly as a nuisance. An interesting concomitant out- 
come of this calibration program is that weak lensing 
surveys will provide constraints on halo structure that 
complement constraints obtained via alternative means 
(e.g., Ref. [imilllllS,!!!,!!!,!!!,!!^). This byprod- 
uct of dark energy studies extends the scientiflc reach of 
future surveys as it may prove useful to inform models of 
the formation of galaxies and galaxy clusters. 

Alongside this success, we emphasize that our results 
must be considered with some circumspection. Though 
it appears to be the dominant effect, we do not expect 
that the influence of baryons is strictly limited to the in- 
ternal structures of halos. In our study we also neglect all 
systematics other than baryonic physics, such as intrinsic 
alignments i37|, iSj, [H, . We consider this exercise to 
be a proof of concept indicating that self calibration is 
achievable at minimal cost and that weak lensing tomog- 
raphy may constrain halo profiles at interesting levels. 

There are several papers in the literature that are 
closely related to our work. We have already mentioned 
that l30th White [1^ and Zhan and Knox ^] studied 
the effects of baryons on convergence spectra using ana- 
lytic models and indicated that they will likely be non- 
negligible, at least over some range of scales. Though pre- 
scient, both of these studies are based on phenomenolog- 
ical models that do not treat all of the effects of baryons 



self consistently, while the study of Zhan and Knox [2^ 
focuses on the influence of hot baryons and neglect bary- 
onic cooling. The result of these studies is that bary- 
onic processes influence lensing observables on scales near 
I ^ 10"^ at a level below that reported in Refs. [27i, i28{. 
Furthermore, it is not entirely surprising that the infor- 
mation within forthcoming lensing surveys suffices to cal- 
ibrate additional physics, such as the influence of bary- 
onic processes. Huterer et al. [l^l studied self calibration 
of multiplicative and additive observational systematics 
and found that these could be calibrated successfully at 
a relatively low cost. Self calibration of systematics, in- 
cluding intrinsic alignment effects is also discussed briefly 
in the appendices of the report of the Dark Energy Task 
Force (DETF, Ref. [4l!|) and such self calibration is in- 
cluded in their analyses. 

We describe our methods in the following section. We 
begin with a review of weak lensing tomography. Next, 
we summarize the salient features of the halo model for 
nonlinear clustering. We then review the Fisher matrix 
formalism for assessing the constraining power of forth- 
coming experiments and mention some of the specifics 
of our implementation, including our fiducial model and 
priors. We close the second section with a variety of 
models that we use to incorporate parameterized bary- 
onic physics into our calculations. 

In Section llllj we present the results of our study. 
We open with a review of forecasts for dark energy con- 
straints in the absence of any unknown baryonic physics. 
The results of this exercise represent constraints in the 
limit of perfect knowledge of halo profiles and correspond 
directly to other forecasts in the literature. This limit 
serves as a useful reference for interpreting subsequent 
parameter constraints. Following this, we give forecasts 
for dark energy parameter constraints after including and 
marginalizing over the contribution of unknown baryonic 
physics. We continue with constraints on halo profiles 
themselves after marginalizing over dark energy. Our fi- 
nal results show dark energy constraints behave as a func- 
tion of external, independent constraints on halo profiles 
that can be included in the analysis. We discuss our re- 
sults, including caveats and avenues for future study, and 
draw conclusions in Section [TVl 



II. METHODS 

A. Weak Lensing Observables 

In this study, we consider utilizing information in weak 
lensing observables only. We model our approach after 
that of Ma et al. [13]. In particular, we consider con- 
vergence power spectra from A^tom tomographic redshift 
bins defined by the photometric redshifts of source galax- 
ies. The A^tom(A^tom + l)/2 observables given by the 
number density-weighted convergence spectra and cross 
spectra for the tomographic bins as a function of multi- 
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pole f , are 



B. A Halo Model Primer 



(1) 



In Eq. H{z) is the Hubble expansion rate, Z?a is 
the angular diameter distance, and P{k, z) is the three- 
dimensional total matter power spectrum at wavenum- 
ber k and redshift z. The lensing weight functions VFi(z) 
specify the tomographic bins. Defining the true redshift 
distribution of source galaxies in the ith photometric red- 
shift bin as dni/dz, the lensing weight functions are given 
by 



= -nMHl{l + z)DA{z) 



Da{z,z') dni 
Da{z') dz'- ^ ' 



The present Hubble rate is i/o and Da{z, z') denotes the 
angular diameter distance between redshifts z and z' . 
Note that the source distribution is not normalized to 
integrate to unity and hence P^l represents the power 
spectra weighted by the product of angular number den- 
sities in the ith and jth photometric redshift bins. 

It is natural to express the photometric redshift bins 
in terms of the total true redshift distribution of source 
galaxies An/dz, and the probability of yielding a pho- 
tometric redshift Zp given a true redshift z, which we 
denote P[zp\z). In this case, the true redshift distribu- 
tion of sources in the zth photometric redshift bin can be 
obtained by integrating between the photometric redshift 
limits of the bin 



dni(z) 
dz 



high 



dz. 



dn(z) 



P(zp|z) 



(3) 



For the purpose of this study, we take the true redshift 
distribution to be 



To evaluate the relations in Section III Al it is neces- 
sary to compute the nonlinear matter power spectrum. 
Though other methods exist and have proven quite suc- 
cessful [H, 113] , we use the phenomenological halo model 
to accomplish this. Elements of the halo model have been 
developed in a number of studies [i^, IH, li?!. [isl . A 
comprehensive review can be found in Re f. [5(]|| . We im- 
plement the halo model as in Rudd et al. [28'| and review 
the salient features below. 

The halo model is predicated on the assumption that 
all matter resides within dark matter halos. The matter 
power spectrum is given by a sum of two terms. 



P(fc) = PlH(fc) + P2H(fc). 



(5) 



The first term is from matter elements that reside within 
a common dark matter halo, while the second term is due 
to elements that reside in distinct dark matter halos. The 
utility of this decomposition in the present context is ap- 
parent. Baryonic processes have a significant influence on 
the structures of individual dark matter halos, but have 
little effect on the large-scale clustering of halos them- 
selves and these ingredients are treated independently in 
the halo model. 

The one-halo contribution is 



Pm{k) 



1 

Pm 



dm m? - — }?{k,m), 
dm 
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where pM is the mean matter density of the universe, m 
is halo mass, dn/dm is the mass function of dark mat- 
ter halos, and \{k, m) is the Fourier transform of the 
halo matter density profile, which we have assumed to be 
spherically symmetric and normalized such that the in- 
tegral of the matter density profile p{r, m) , over all space 
is unity, 47r / dr r'^p{r, m) = 1. The two-halo term is 



dn(z) 
dz 



2nz§ 



exp[-(z/zo)^] 



(4) 



with zq ~ 0.92, which fixes the median survey redshift to 
^mcd — 1- The parameter n represents the total density 
of source galaxies per unit of solid angle. For simplicity 
and concreteness, we model the photometric redshift dis- 
tribution as in Model I of Ref. [i^. To be specific, we 
take P{zp\z) to be a Gaussian centered at Zp = z with dis- 
persion az = 0.05(1 -I- z). We define tomographic bins to 
be evenly spaced in redshift from z = to z = 3, putting 
whatever fraction (of order ~ 10^'*) of galaxies that lie 
beyond z = 3 into the highest redshift tomographic bin. 
We number our tomographic bins sequentially so that 1 
designates the lowest redshift bin and TVtom labels the 
highest redshift bin. Unless we state otherwise, we re- 
port results with A^tom = 5 as we find little relative 
improvement in parameter constraints with finer tomo- 
graphic binning, in agreement with previous work [42] . 
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dfi 

dm m- — X{k,m)bh{m) 
dm 



(7) 

where P'™(fc) is the matter power spectrum given by lin- 
ear perturbation theory and b\^{m) is the mass-dependent 
halo bias. We calculate the linear matter power spectrum 
using the fitting formulae of Eisenstein and Hu ISlll with 
the appropriate modifications for dark energy [52] . The 
halo model has known shortcomings (e.g., [53|), but we 
model the matter power spectrum in this way because 
it isolates the quantities we aim to study, namely the 
mean density profiles of halos of mass m, p(r, m) and 
their Fourier transforms X{k,m). 

Motivated by the findings of Ref. [11], we model the 
average mass density of halos using the density profile of 
Ref. El, 



p{r) oc 



1 



(cr/i?200m)(l + Cr/i?200n 



(8) 
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truncated for r > i?200m- The parameter c describes 
the relative concentration of mass toward the halo cen- 
ter. i?2oom gives the extent of halo in the sense that the 
mass within i?200m is the halo mass m. We define ha- 
los as spherical objects within which the mean density 
is 200pM, so that the mass and radius are related by 
m = 47r(200pM)^200m/3. This convention is chosen to 
be consistent with the halo bias prescription that we use 
(see below) . The density profile of Eq. ([8|) has a Fourier 
transform [4811 



A(fc,c) 



sin(r;)[Si([l-fc]r;)-Si(r7)] 



fic) ^ 

-fcos(77)[Ci([l + c]7?)-Ci(?7)] 
sin(77) I 



CJ77 



(9) 



where f{x) — ln(l + x) — x/{l + x), r] = fci?200m/c, and 
Si(a;) and Ci(a;) are the sine and cosine integrals respec- 
tively. 

In the case of dissipationless physics, we set the mean 
halo concentration as a function of mass and redshift ac- 
cording to a power law [11] (see also Refs. [H, 113, HI] and 
in particular Ref . [m^] for a recent study of concentrations 
in the Millennium Simulation [85]) 



c^^"(m,z) = ll[m/m,,o]""-'(l + ^)" 



(10) 



where to*,o is the mass of a typical object collapsing at 
z = {nii^fi ~ 2.2 X 10^^ /i~^Mq in our fiducial cosmol- 
ogy, see below). The slight difference between Eq. ((TU]) 
and the relation in Ref. [55| is due to the conversion be- 
tween the the halo mass definition used in that study and 
the one we use Throughout this study we ignore 

the distribution in concentrations at fixed mass which is 
approximately log- normal [s^ . For concentration distri- 
butions computed from TV-body and hydrodynamic sim- 
ulations, the relative influence of the spread in concentra- 
tions on convergence s pec tra is completely negligible on 
the scales we consider [6ll . l62j . For simplicity, we model 
the halo mass function and the halo bias using the rela- 
tions of Ref. [63| . We reiterate here that when we refer to 
halo density profiles and halo concentration parameters, 
we mean the effective density profiles and concentrations 
of the composite systems of dark matter and baryons. 



C. Fisher Matrix Analysis, Experimental 
Specifications, and Cosmological Parameters 

We quantify the constraining power of observables us- 
ing the Fisher information matrix. The Fisher matrix 
components are a sum over observables, 

= x;(2^+i)/.k..E?^[^"']AB^+^ir' (11) 



where the Oa represent the iVTOM(-^TOM + l)/2 observ- 
ables of Eq. Il]) and the index A runs over these observ- 
ables. Cab are the components of the covariance matrix 
of observables, and the pi represent the parameters of 
the model. The sum over multipoles runs from some 
minimum multipole £min fixed by the sky coverage of the 



2/s 



-1/2 



where 



experiment. For simplicity, we set ^uim — ^jg^y 
/sky is the fraction of the sky covered by the experiment; 
however, this choice is unimportant because constraints 
are dominated by considerably higher multipoles in all 
cases of interest. The multipole sum runs to some max- 
imum multipole £max, which is typically fixed between 
^max ~ 1000 — 3000 so as to remain in the regime where 
several assumptions, including that of Gaussian statis- 
tics, are valid Bl], [gj, [H, [H, [63] • Furthermore, our mod- 
eling of baryonic effects is certainly not valid at signif- 
icantly higher multipoles, so it is sensible to keep ^max 
in this range for our purposes. We elaborate on this in 
section III Dl 

The inverse of the Fisher matrix approximates the 
parameter covariance locally about the maximum like- 
lihood. As such, the measurement error of the ith pa- 
rameter is cr(pi) = [F^^]ii. The second term in Eq. (fTTI) 
incorporates any Gaussian priors on the model parame- 
ters. In most of what follows, we assume modest priors on 
each parameter individually so that = Sij/af, where 
6ij is the Kronecker S symbol. We discuss our parameters 
and priors in more detail below. 

The Fisher matrix formalism also provides an estimate 
of parameter biases due to unknown systematic offsets 
in observables. Let O^^^ be the difference between the 
true observable and the perturbed observable due to some 
systematic effect. In the limit of small systematic offsets, 
the parameters extracted from this set of observables will 
be biased by 



9pj 
(12) 

Observed spectra contain both a term due to signal 
and a noise term, 



i^(£) = f«(^)+ni5ij(7^). 



(13) 



where (7^) is the intrinsic source galaxy shape noise, and 
rii is the surface density of sources in the ith tomographic 
bin and can be obtained by integrating Eq. ([3]) over all 
redshifts. The covariance between observables P^' and 

pkl 



IS 



(14) 



where the i and j map to the observable index A, k and 
1 map to B, and we have assumed the stati stics of the 
convergence field to be Gaussian. We fix (7^) = 0.2 
and study constraints from several forthcoming surveys 
defined by the parameters /sky and n, the total effec- 
tive surface density of source galaxies on the sky. For 
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specificity, we study the limits expected from an experi- 
ment like the Dark Energy Sm'vey (DES)|86i] with /sky = 
0.12 and n — 15/arcmin^, a future space-based mission 
like the SuperNova/ Acceleration Probe (SNAP) [87] with 
/sky = 0.025 and n = 100/arcmin^, and a future ground- 
based survey like the Large Synoptic Survey Telescope 
(LSST)^'88l| with /^ky = 0.5 and n = 50/arcmin^ 

We consider cosmologies defined by seven parame- 
ters. Three of these parameters describe the dark en- 
ergy. The other four parameters and the values they 
take in our fiducial cosmological model are the total mat- 
ter density ujm = ^uh^ — 0.128, the baryon density 
lub = 0.0223, the amplitude of curvature fluctuations at 
k = 0.05 Mpc^^ = 2.04 x 10"^ (we actually vary 
In A|j), and the power-law index of the primordial power 
spectrum rig = 0.958. This fiducial model is motivated 
by the Wilkinson Microwave Anisotropy Probe three- 
year results 0. We take very modest priors on these 
four parameters, (t(wm) = 0.007, (j{ujb) = 1-2 x 10~^, 
(T(lnA^) = 0.1, and a{ns) — 0.05, which are compara- 
ble to contemporary constraints (e.g., Refs. [1, We 
describe the dark energy by its present energy density pa- 
rameter r^DE = 0.76, and an equation of state parameter 
that varies with scale factor as w{a) = wg + {1 — a)?i'a. In 
our fiducial model, wq = —1 and = and we assume 
no prior information on any dark energy parameter. As a 
reference, our fiducial model implies that the root-mean- 
square matter density fluctuations on a scale of 8 /i"^Mpc 
is as — 0.76. 

In addition to the marginalized errors in the cosmolog- 
ical parameters, we also quantify the constraining power 
of surveys by the figure of merit advocated in the report 
of the Dark Energy Task Force pl| (see Ref. for a 
similar suggestion). The DETF figure of merit is the in- 
verse of the area of the marginalized 95% ellipse in the 
WarWo plane divided by tt. Letting a and h denote the 
lengths of the principal axes of this ellipse, the figure of 
merit is = 1 /ab. Many studies also quote a pivot scale 
factor Cp, and a pivot value of the dark energy equation 
of state parameter w-p = w{ap) [ill, HI, HI]. The pivot 
scale factor is the scale factor at which the uncertainty 
in w{a) is minimized, and is given by 



1 



\F- 



The pivot equation of state parameter is simply 



Wo + (1 - ap)wa. 



(15) 



(16) 



Taking Wa. along with as the parameters describing 
the dark energy equation of state is convenient because 
the principal axes of the error ellipse lie along these co- 
ordinates. Moreover, the transformation is linear and 
preserves the area of the ellipse. In terms of the compo- 
nents of the Fisher matrix [Eq. (jlip ]. the variance in Wp 
is given by 



a\wp) = [F-% 



-112 



(17) 



The 95% confidence level in two dimensions occurs at 
~ 2.48fT, so that the area of the 95% ellipse is ~ 
6.177r(7(wa)(T(wp). The DETF also quote values of 
[cr(i(;a)c(wp)]~^ in their summary tables 4l|. Clearly, 
our figure of merit T , is related to the numerical values 
quoted in the tables of the DETF report [cr{w-p)a{'Ws.)\^^ , 
by 



0.162 



(18) 



We intend to focus on the relative scaling of the figure of 
merit, so the details of the constants of proportionality 
here are of minimal importance. 



D. Modeling the Effects of Baryons with Modified 
Halo Concentrations 



Rudd et al. [28| demonstrated that the influence of 
baryons on the matter power spectrum is primarily due 
to the modified structure of individual dark matter halos. 
In fact, the modification is simple. Rudd et al. [28j found 
that the halos in their hydrodynamic simulations followed 
a concentration-mass relation that was significantly dif- 
ferent from the relation derived from iV-body simulations 
[e.g., Eq. pH]) ]. Rudd et al. [28'] found that the descrip- 
tion of halo structure as NEW halos with boosted con- 
centrations was valid for halo-centric radii greater than 
approximately 0.04i?2oom [111- Upon accounting for this 
difference, these authors found that they could model 
convergence power spectra derived from their simulations 
accurately out to multipoles as high a,s i ^ 6000. How- 
ever, as noted in both Ref. [23| and Ref. [2^, hydro- 
dynamic simulations are not yet up to the task of pro- 
viding accurate predictions for the properties of galax- 
ies. So, while the net effect of galaxy formation on the 
power spectra may be well described by a modified, ef- 
fective c(m, z) relation, this effective c(m, z) relation is 
not known. Consequently, it seems natural to calibrate 
the effective concentration-mass relation within the ob- 
servational data itself and we explore several parameteri- 
zations of varying complexity for this relation. In the fol- 
lowing section, we present results where both cosmology 
and the parameters of the c(m, z) relation are determined 
simultaneously from observational data. 

Neglecting baryonic physics, the relation between con- 
centration and mass is known to be well characterized by 
a power law [Eq. pU]) . see Ref. [El, [5^]. Consequently, it 
seems natural to study power-law relations of the form 



c(to, z) = co[TO/m^^o]"(l -I- z) 



(19) 



Further, it is reasonable to suppose that the redshift de- 
pendence is not modified from the standard A^-body re- 
lation [Eq. pI7|) ] because the majority of baryonic con- 
densation is thought to have occurred at 1 and be- 
cause the redshift scaling in the standard case arises sim- 
ply because virial radii grow as -R200m oc (1 -I- z)~^ in 
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the absence of accretion. This scahng arises because ha- 
los are typicaUy defined as an overdensity with respect 
to the mean density and the mean density dilutes as 
Pm (1 + while the structure of the bound object 
remains unchanged. We expect this second feature to 
be unmodified when including baryonic physics. Nev- 
ertheless, we study concentration relations of this form 
both with the redshift dependence fixed to /3 = 1 and 
with /3 free. We choose a fiducial concentration rela- 
tion that is enhanced so that Cq = 15 and a — —0.1. 
This represents a 36% enhancement in halo concentra- 
tions while Rudd et al. [2^ report a slightly larger en- 
hancement in concentrations in their DMG_SF simulation 
for halos near m a few x 10^^ h~-^MQ. As we demon- 
strate in Section imi it is the halos in this mass range that 
most infiuence convergence power spectra. We choose 
an enhancement slightly below the most relevant result 
from Ref. [28| because the simulation results likely over- 
estimate the effective concentrations due to the well- 
documented overcooling problem of cont emp orary hydro- 
dynamic simulations (ZO, llll, [zl, [zl, [zi S IH ■ 

Though concentration relations normalized at m — 
m*jO and redshift z = are conventional in the litera- 
ture on halo structure, this is likely not to be the most 
appropriate representation in the present context. Halos 
near m ~ 10'^'' /i~^Mq make the largest contribution to 
convergence power spectra [1^, |z3 , and lensing weights 
are considerable at redshifts between 0.2 < z < 1, so 
it is natural to suppose that lensing will be most sensi- 
tive to the structure of ~ 10^'* h~-^MQ halos in a broad 
range of redshifts. In analogy with the the pivot scale 
factor and pivot equation of state parameter described 
in Section fH Cl for the dark energy, we may define a pivot 
redshift Zpiv, and pivot mass rripiv, for the c(m, z) rela- 
tion. We set the Zpiv and the redshift and mass 
at which the fractional uncertainty a{c)/c is minimized. 
The pivot concentration is Cpiv = c(mpiv, Zp 
tional uncertainty at the pivot concentration is 



<-piv)- The frac- 



a^cpw) _ [F- 



-I- ln^(mpiv/m*,o)[^ ^]ac 



ln2(l + Zpiv)[F-i]^/3 

21n(TOpiv/m*,o) ri;,_ii 



21n(l + Zpiv), 1, 

-[^ Jco/3 



Co 



- 21n(mpiv/m,^o)ln(l + Zpi^)[F-\p.{2Q) 
The pivot redshift is 

1 [F~^]aa[F~'^]col3 - [F~%oa[F~'^]a0 



ln(l -I- Zpiv) 



Co [F-^U[F-^00-[F-m 



(21) 

and the pivot mass is 

1 [F-^00[F-^]coc. - [F-\MF-'U 



ln(mpiv/m*^o) 



CO [F-m,-[F-^^^[F~^^,p 



(22) 



In the following sections, we will quote pivot masses and 
redshifts and errors in Cpiv. 

The c(m, z) relation of Eq. (jl9p may be general enough 
to encompass most plausible behavior, particularly be- 
cause the convergence spectra are most sensitive to a 
narrow range of halo masses near m ^ lO^"' /i^^Mq for 
£ ^ 10^ [13, [z3 (also see Sec. HH. However, it is not 
unreasonable to suspect more complicated behavior. For 
example, the mass-to-light ratios of galaxy, group, and 
cluster halos are known to vary in a non-monotonic fash- 
ion with total mass (e.g., Ref. [78, 79, 80] ) with a mini- 
mum near m ~ 10^^'^ /i^^Mq. As this minimum reflects 
some maximum efficiency for galaxy formation, it is nat- 
ural to suspect that the relative boost in halo concentra- 
tions peaks at this halo mass and declines on either side 
of it. To hedge against the possibility of more compli- 
cated behavior than that in Eq. ([19]), we also explore a 
concentration relation with significantly more freedom. 

In this more general relation we define values of concen- 
tration at Nc values of log(TO/ ^."^M©) evenly spaced in 
the log of the halo mass between 11 < log(TO/ h~^MQ) < 
16. We have verified that extending this range to lower or 
higher masses has little effect as these masses contribute 
too little to the convergence spectra on scales of interest 
[26l Itt} . The concentration values specify the mean halo 
concentration at the center of each mass bin and we spline 
interpolate between these values to enforce smoothness in 
the relation. This choice refiects our prejudice that the 
average concentration should always be a very smooth 
function of mass. In addition, we expect that at very high 
masses and very low masses baryonic condensation and 
galaxy formation should be inefficient, so we require that 
the standard relation [Eq. (|10p ] be recovered for masses 
log(m/ /i-iM©) < 9 and log(m/ /i-^Mq) > 18. We en- 
force this by matching the spline onto the standard rela- 
tion at these endpoints. In practice, the choice of bound- 
ary condition for the c(to, z) at low and high masses is 
unimportant. We refer to our parameter set by subscript- 
ing a "c" with the logarithm of the mass at which this 
concentration is specified. For example, in a model with 
Nc — 5, our concentration parameters are C11.5, C12.5, 
C13.5, C14.5, and C15.5. We implement a variable redshift 
dependence as c(m, z) = c{m, z = 0)/(l + z)^ . 

To give the reader a sense of the influence of these 
changes on observables, we show the relative differences 
in the observable convergence spectra in Figure [TJ For 
concreteness, we plot spectra in the power law c(to, z) 
model with co = 15 compared to the spectra computed 
using the standard relation of Eq. ([TUj) derived from dis- 
sipationless simulations. We have used a model with five 
tomographic redshift bins, but for the sake of clarity we 
plot only three observables. These correspond to the 
spectra of the second and third tomographic bins and 
the cross spectrum between these bins. As a reminder, 
the photometric redshift limits of the second and third 
bins are 0.6 < Zp < 1.2 and 1.2 < Zp < 1.8 respectively. 
As one would expect, any effect on concentrations be- 
comes quite large at 10'^. However, it is worth noting 
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FIG. 1: The differences in convergence spectra for the fiducial 
model of boosted power-law concentrations {thin lines) rela- 
tive to the standard case with concentrations derived from 
dissipationless simulations (thick lines). In the top panel we 
show the power spectra explicitly and in the bottom panel 
we show the change in observables scaled to the standard 
TV-body concentration case. We use five tomographic bins 
throughout most of this work and we have dropped the num- 
ber weighting of spectra for the purposes of this demonstra- 
tion only. For simplicity we show only the behavior of three 
observables, the spectra in the second and third redshift bins 
P^^/n^ and Pj^'^/n^, and the cross spectrum P^^/n2n3. The 
behavior of other observables is qualitatively similar. The 
shaded bands show, from outermost to the innermost at left, 
the variance in P^^ for SNAP, DES, and LSST in bands of 
width A£/i = 1/10. Note that the sensitivity of SNAP over- 
takes that of DES by £ ~ 600 as SNAP is far deeper than 
DES. We show this by extending the upper limit of the shaded 
SNAP band as a solid line below the filled DES region. 



that the effect is not negligible even at much lower mul- 
tipoles and that all of the experiments we consider would 
be sensitive to such effects even at multipoles as low as 
several hundred. 



III. RESULTS 
A. Dark Energy Constraints and Parameter Biases 

We begin our presentation of results with parameter 
constraints in the standard case of no modifications to 
halo structure from baryonic processes. This case corre- 
sponds to perfect knowledge of halo density profiles and 
has been considered numerous times in the literature al- 
ready. We include this partly to give a frame of reference 
for the results that follow. More importantly, we show 
dark energy parameter constraints as a function of the 
maximum multipole fmax used in the sum in Eq. (jlip 
so that one can compare the relative merits of marginal- 
izing over an unknown halo concentration-mass relation 
with simply ignoring the high multipole moments where 
baryonic effects are most severe. Throughout this sec- 
tion, we label the uncertainties in parameter p obtained 
by including multipoles out to imax = £ as oi(p). 

Constraints on dark energy parameters in the ab- 
sence of any unknown baryonic processes as a function 
of ^max are shown in Figure [U where we have taken 
.^TOM = 5 as prescribed by Ref. [i^. Unless other- 
wise noted, this choice pertains to all results on pa- 
rameter constraints. With modest prior information, 
weak lensing tomography alone can constrain dark en- 
ergy parameters very well. Using information out to 
an £niax = 3000, the marginalized constraints from 
DES, SNAP, and LSST respectively are crfj'oo (^^de) = 
0.023,0.018,0.0066, ctIto^Iw^o) = 0.26,0.19,0.069, and 
c^3(X)o(w'a) = 0.74,0.56,0.21. Constraints from LSST are 
the most stringent because the experiment combines ex- 
tensive sky coverage with depth. SNAP and DES are 
comparable (considering only weak lensing information) 
over this range of scales. SNAP, being extraordinar- 
ily deep, can exploit information in very high multi- 
poles while the constraints from DES overtake those from 
SNAP if only multipoles 400 are used because of its 
larger sky coverage. Hereafter, we plot our results in- 
cluding self-calibration of baryonic processes relative to 
these standard results which represent the limit of perfect 
knowledge of internal halo structure. 

Both Jing et al. ^ and Rudd et al. [Ill pointed out 
that the modifications of the convergence power spec- 
tra due to baryonic physics are likely to be large and to 
extend over a wide range of multipoles, but neither of 
these studies translated this systematic into its influence 
on cosmological parameter constraints. In particular, it 
is interesting to examine the bias in dark energy parame- 
ters if baryonic physics are neglected. If the bias is small 
compared to the statistical uncertainties in Fig. ^ then 
it need not be a major concern because it is unlikely to 
lead to rejection of the true cosmology based on these 
data. 

We show the biases in the inferred values of wq and 
as a function of ^max in units of their uncertainties 
in Figure [31 To compute the biases shown in Fig. [31 
we assumed the true power spectra to be given by the 
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FIG. 2: Standard marginalized dark energy parameter con- 
straints as a function of maximum multipole ^max. The left 
panels show the la uncertainties in the three dark energy pa- 
rameters Qde, Wo, and Wa from the top down. The bottom 
panel shows the dark energy figure of merit .7^ as a function of 
maximum multipole. The solid line corresponds to an LSST- 
like experiment, the dashed line corresponds to a SNAP-like 
experiment, and the dotted line corresponds to an experiment 
similar to DES. The right panels show the same quantities 
scaled by the uncertainties obtained at ^max = 3000 in order 
to show the relative scaling of dark energy constraints with 
^max more clearly. 



fiducial power-law concentration model with cq — 15, 
which is a ~ 36% enhancement relative to the standard 
concentration-mass relation, and employed the approxi- 
mation of Eq. (|12p . The magnitudes of the biases show 
sharp minima in many cases. These represent changes 
in the signs of the biases. Generally, 5{wo) > and 
S{wa) < at the lowest multipoles. Note that in all 
cases biases become large compared to the uncertainties 
in the parameters, so the approximation used to compute 
them breaks down. These large biases are unlikely to be 
the product of any analysis. Rather the team performing 
the analysis would notice something amiss because the 
best-fitting models would be relatively poor fits to the 
data and this would be evident in some criterion used to 
assess fit quality, such as the P^r degree of freedom. 

The details of the parameter biases are a strong func- 
tion of the choice of fiducial cosmology. This is because 
both the signal-to-noise ratios and the scales at which 
nonlinear effects become important are strong functions 



of the fiducial model. We illustrate this in Fig. [3] by 
showing biases computed about the fiducial cosmologi- 
cal model of Ma et al. [i^l , which more closely resembles 
the cosmology from the Wilkinson Microwave Anisotropy 
Probe first year results ^81,]. The most important differ- 
ences between our fiducial model and the fiducial model 
of Ma et al. ^ are that Ma et al. [i^l took = 1 
(compared to Us = 0.958) and a slightly higher power 
spectrum normalization, implying a larger mass variance 
on 8 h-^Mpc scales of as ^ 0.9. In the Ma et al. ^ 
fiducial model, parameter inferences are more sensitive 
to small systematics on quasi-linear and nonlinear scales 
for two reasons. First, the signal is relatively higher. Sec- 
ond, nonlinear effects become important on larger scales 
(smaller £) in models with more power. Both of these 
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effects are reflected in Fig. |31 Near the Ma et al 
fiducial model with relatively greater power, biases are 
typically larger and they become important at smaller 
wavenumbers. Throughout the remainder of this section, 
we give our primary results derived about the fiducial 
model of Section fll C\ however, we often compare these 
to results derived about the fiducial model of Ma et al. 
42] in order to indicate the importance of the choice of 
fiducial model. 

As the details of the fiducial model have a significant 
influence on parameter biases and the linearized relation 
of Eq. (fT^ should break down for biases large compared 
to their uncertainties, it is hard to make concrete state- 
ments about parameter biases. Nevertheless, Fig. [3] is 
an indicator that such biases may be signiflcant, even 
at multipoles £ < 10"^. LSST will provide the tightest 
constraints on dark energy parameters over this range 
of scales, primarily because it is extremely wide, yet 
this makes LSST most susceptible to unknown baryonic 
physics in a relative sense. In the absence of reliable 
and robust predictions for power spectra in the presence 
of baryons or any other method for dealing with the in- 
fluence of baryonic processes, one could safely force the 
biases to be relatively small compared to the uncertain- 
ties by disregarding multipoles above some imax- This 
and more elaborate strategies for removing some of the 
small-scale information in weak lensing tomography have 
been studied by Huterer and White [84I . Fig. [3] suggests 
that this would require setting ^max to no greater than 
many hundreds, though the details depend upon both the 
experiment (they are most stringent for LSST) and the 
underlying cosmology. Comparing with Fig. [2l this cor- 
responds to a degradation of dark energy parameter con- 
straints by a factor of ^ 2 — 3 or more relative to the con- 
straints achieved by setting ^max = 3000 in the standard 
case of perfect knowledge of halo density profiles. Corre- 
spondingly, the decrease in the DETF figure of merit is 
approximately a factor of ~ 4 — 7. We note that these ef- 
fects are more important in models with more small-scale 
power, so that about a cosmology with a higher value of 
(Tg ~ 0.9, dark energy parameter constraints can be de- 
graded by more than a factor of three and the figure of 
merit can be decreased by nearly a factor of ten. In what 
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FIG. 3: Biases in inferred dark energy parameters due the 
ignorance of baryonic physics. We show biases by assuming 
the true model to have a power-law concentration-mass rela- 
tion with Co ~ 15 instead of co = 11 as in the case of models 
that neglect baryonic physics. Results for DES, SNAP, and 
LSST are shown as labeled in the upper left portion of the 
figure. In all cases, we show bias relative to the la measure- 
ment uncertainty in the relevant parameter. The top panels 
show the biases for wo and the bottom panels show the biases 
for Wa,- Notice that we show biases about both the fiducial 
model that we use in this work {left column) and the model of 
Ma at al. [i^ which has significantly more small-scale power 
{as — 0.9, right column). The minima in the magnitudes of 
the biases reflect changes in sign. Generally S{wo) > at the 
lowest multipoles and 5{wa) < at the lowest multipoles. In 
the case of LSST, the first change in sign in S{wo) occurs at 
^max < 300, so that 5{wo) < at ^max = 300. The approxima- 
tions used to compute the parameter bias [see Eq. (|12p ] break 
down in all cases as the biases become larger than their un- 
certainties at high multipoles. These estimates indicate that 
the biases are likely to be non-negligible. 
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FIG. 4: The relative degradation in parameter constraints 
as a function of maximum multipole used to derive con- 
straints for the case of a power-law concentration mass re- 
lation [Eq. (|19p ]. The vertical axes are the parameter con- 
straints in the new concentration model divided by the con- 
straints in the standard case. In the left panel, we show 
constraints in a model with fixed redshift dependence. The 
thick lines show constraints about a fiducial model with con- 
centrations enhanced relative to the standard A''-body case 
(co = 15) while the thin lines show constraints about the A''- 
body (^^{m^z) (cq = 11) as fiducial model [Eq. US}]. In 
the right panel we show constraints in a model in which the 
power-law index describing the redshift dependence of con- 
centrations is variable. The fiducial model in this case and in 
all subsequent plots has enhanced halo concentrations rela- 
tive to the dissipationless case. Note that the dynamic ranges 
on the vertical axes are different in the left and right panels. 
In all panels, the solid line represents LSST, the dashed line 
represents SNAP, and the dotted line represents DES. 



follows, we show that internal self calibration of baryonic 
effects may be a viable alternative to this degradation in 
parameter constraints. 

B. Dark Energy Constraint Degradation: 
Power-Law Halo Mass-Concentration Relations 

We consider the degradation in parameter constraints 
sequentially from the power-law c(m, z) model to the 
more general cases in order to give a sense of the grad- 
ual degradation from the additional degrees of freedom. 
The relative degradation in dark energy parameters in 



the power-law c(m, z) models is shown in Figure [H The 
left panel of Fig. H] shows constraints with a fixed red- 
shift dependence and the right panel shows constraints 
with a variable redshift dependence parameter /3. This 
figure shows that the degradation of dark energy param- 
eter constraints is relatively small and that it behooves 
one to model the unknown baryonic effects and continue 
to use high multipole moments rather than cut out infor- 
mation above some relatively low £max- Provided there 
is some compelling justification for these relatively re- 
stricted concentration models, we can eliminate the bias 
caused by unknown halo structure while dark energy con- 
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straints will be degraded by only '--^ 15 — 20% or less in 
all cases. This cost is relatively low compared to the 
cost of using only multipoles below several hundred to 
constrain dark energy and, as we discuss below, the nui- 
sance parameters of this analysis are physical parameters 
of interest in their own right. 

Aside from this overriding point, there are several other 
features worthy of note. First, notice that there are 
instances (depending upon multipole range and exper- 
iment) in which the model with extra freedom provides 
better constraints than the standard model with known 
halo concentrations. This is due to our assumption that 
concentrations are enhanced (we take cq = 15 in our 
model with baryons compared to cq = 11 as given by 
A^-body simulations) resulting in increased signal over a 
wide range of scales (See Fig. [1]). The thin lines in the 
left panel of Fig. [4] show constraints about the standard 
c(m, z) relation with cq = 11 as fiducial model. In this 
case, there is no boost in signal and the constraints are al- 
ways poorer than in the model with perfect knowledge of 
c(m, z), though not by much. Hereafter, we quote results 
for a fiducial model with enhanced halo concentrations 
relative to the standard iV-body values. In all cases, the 
choice of fiducial cosmological parameters affects param- 
eter constraints at levels similar to those shown in the 
left panel of Fig. g) 

In addition, observe that in the panels of Fig. [?] the 
relative degradation has a peak at multipoles of a few 
thousand with some variation depending upon the pa- 
rameter of interest and the experiment. The peaks are 
more prominent in the right panels. In fact these fea- 
tures would also be apparent in the left panels of the 
figure, but are mostly omitted as they lie beyond our 
choice of £max — 3000. The origins of these features are 
simple to understand. The effect of concentrations is 
predominantly a high-^ phenomenon. As one approaches 
i ~ 2000 from the low multipole side, the unknown con- 
centrations are degrading dark energy constraints rapidly 
because their effects are just beginning to appear and can 
at some level be mimicked by the dark energy param- 
eters. However, as one continues to higher multipoles 
(£^ 3 X 10^), the concentration effects become increas- 
ingly distinguishable from those attributable to dark en- 
ergy parameters, these degeneracies are broken and con- 
straints improve to higher multipoles thereafter. The 
low-^ feature in these plots is due to a similar degen- 
eracy primarily between the power-law index a and the 
dark energy equation of state parameters. Such features 
caused by the gradual elimination of degeneracies with 
the inclusion of information from ever higher multipoles 
are generic in the concentration models that we explore 
and are present in subsequent plots as well. In all cases, 
these "rolling hills" represent the sequential introduction 
and elimination of parameter degeneracies and we do not 
expound upon the details of each degenerate direction in 
this study in the interest of brevity. 

It is important to stress that tomographic information 
is a crucial ingredient in minimizing parameter constraint 
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FIG. 5: The degradation in dark energy parameter con- 
straints without tomographic information. This figure shows 
the relative errors compared to the case of perfect knowledge 
of effective halo concentrations as a function of ^max as in 
Fig. 131 but with no tomography. To illustrate this point, we 
use the power law c(m, z) model with variable redshift index 
l3. For the sake of clarity, we show only the case of LSST 
here, but the loss of constraining power is dramatic for all 
experiments. 



degradation in the case of unknown halo concentrations. 
For example, eliminating tomographic information (set- 
ting A^TOM = 1) dramatically reduces the ability of such 
surveys to calibrate for unknown halo concentrations. We 
find that in the case of the power-law class of concen- 
tration models with variable redshift index /3, eliminat- 
ing tomographic information in both the concentration- 
marginalized results and the A^-body standard causes 
the relative degradation to rise from the ~ 15% level 
as shown in Fig. 2] to a factor of ~ 5 for Wa and a fac- 
tor of ~ 10 for SIde and wo- Note that this degradation 
is with respect to an already degraded baseline with no 
tomographic information. Figure [5] shows an example 
of the degradation in dark energy parameter constraints 
without tomography for the LSST experiment. More- 
over, note that this level of dark energy information loss 
extends to multipoles as low as ^max ~ 600. On the other 
hand, we find that constraints on dark energy parameters 
effectively saturate at A^tom = 5 getting only very grad- 
ually better with increasingly fine photometric redshift 
binning. For instance, increasing to A^tom = 10 only 
gives a few percent decrease in dark energy parameter 
errors. 

The importance of tomographic information has both 
positive and negative aspects. On the negative side, it 
stresses the need for accurate photometric redshifts and 
photometric redshift distributions with well-understood 
properties. However, it suggests that the ability of fu- 
ture surveys to self calibrate does not rely solely on the 
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shapes of the power spectra produced by modified bary- 
onic physics being non-degenerate with those arising due 
to dark energy effects. To be more specific, this also sug- 
gests that this success is not limited to very restrictive 
c(m, z) prescriptions. For example, if the power spectra 
shapes obtained by modifying halo concentrations were 
in no way degenerate with the shapes obtained by vary- 
ing dark energy parameters, eliminating tomographic in- 
formation would not lead to further degradation in pa- 
rameter constraints because even in this case the non- 
degenerate parameters could not compensate for each 
other. That constraints do significantly degrade in the 
case of no tomography suggests that dark energy and 
concentrations may lead to degenerate power spectra but 
that tomography helps distinguish the two. Tomographic 
information effectively provides several handles on the 
relative distances to the sources in the different photo- 
metric redshift bins [Eq. ^] and enables one to dis- 
tinguish between contributions to the shear from high 
redshift and low wavenumber from contributions at high 
wavenumber and low redshift, both of which may con- 
tribute at a single multipole. These internal consistency 
checks within weak lensing tomography render the dark 
energy and concentration parameters relatively distin- 
guishable. Already, this suggests that the more general 
concentration relations we describe below will not lead 
to drastic dark energy parameter degradation. 



Constraint Degradation: Binned Halo 
Mass-Concentration Relation 



We now present parameter constraints in the context 
of our more general, binned c(m, z) models. Based on 
the discussion of the previous section, there is cause to 
be optimistic that self calibration with a more general 
relation for halo concentrations will quite successful. In- 
deed, Figure [6] reveals this to be the case. Maximal 
parameter degradation is ~ 30% or less in all cases, so 
again self calibration performs better than the more dras- 
tic measure of disregarding the information contained 
in multipoles larger than a few hundred. Figure [6] de- 
picts the case of a five-bin concentration-mass relation 
for concreteness; however, we have found that the level 
of degradation is relatively robust to the number of bins 
used. In fact, the level of degradation remains below a 
worst case degradation of 50% in f^DE and Wa and 40% 
in Wa for LSST and is below 40% in all parameters for 
DES and SNAP using 12 bins over the same mass range 
(11 < log(m/ H'-^Mq) < 16). Given our prejudice that 
the concentration-mass relation should be a very smooth 
function of halo mass, this suggests that self calibration 
can successfully distinguish the variety of such relations 
realized in the universe while providing stringent con- 
straints on dark energy. 

The features in Fig. [S] are qualitatively similar to those 
in Fig. 21 The fact that there are multiple features in this 
case reflects the additional freedom in the concentration 
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FIG. 6: The relative degradation in parameter constraints as 
a function of maximum multipole used to derive constraints 
for the case of the binned concentration mass relation. The 
vertical axes are the parameter constraints in the new concen- 
tration model divided by the constraints in the standard case. 
In the left panel, we show constraints in a model with fixed 
redshift dependence. In the right panel we show constraints 
in a model in which the power-law index describing concen- 
tration redshift dependence is variable. In all panels, the solid 
line represents LSST, the dashed line represents SNAP, and 
the dotted line represents DES. 



relation. Each concentration bin influences convergence 
spectra over a particular range of scales. The largest ha- 
los contribute on the largest scales, while the smaller ha- 
los become increasingly important with increasing multi- 
pole [H, [z3|. In this case, the concentration parameters 
begin to be constrained gradually, from highest to lowest 
mass, as one increases ^max- Aside from these features, 
the primary aim of Figure[B]is to illustrate the important, 
global point that dark energy parameter constraints are 
robust to very general concentration relations. Being the 
most sensitive instrument over the range of scales we con- 
sider, LSST is most sensitive to these effects, but even in 
this case the cost of avoiding the potential bias shown in 
Fig. [3] is low. 
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D. Constraints on Concentrations from Weak 
Lensing Tomography 

To this point, we have considered the unknown struc- 
ture of halos to be a contaminant that degrades our 
abihty to study the cosmic dark energy. However, ob- 
servational constraints on concentrations are interesting 
in their own right. Indeed, we were driven to consider 
general parameterizations of halo structure because the 
baryonic physics that govern galaxy formation and influ- 
ence the structures of halos is very uncertain. Direct con- 
straints on concentrations may help to inform models of 
galaxy formation. Furthermore, constraints derived from 
a weak lensing survey as described above, would not be 
subject to the same selection as studies of galaxy groups 
and clusters or studies that rely on selecting sample mem- 
bers or stacking systems according to some member prop- 
erty. Therefore, these constraints can complement other 
techniques even if they are not directly competitive with 
other methods. We now focus on the constraints on halo 
concentrations that can be derived from weak lensing to- 
mography. 

Consider first constraints on the effective halo mass- 
concentration relation in the power-law concentration 
mass relation. Figure [7] shows two-dimensional con- 
fidence contours for the concentrations of halos after 
marginalizing over all other parameters, including dark 
energy. For convenience, the pivot masses and redshifts 
are listed in the upper right portion of this plot. Con- 
centrations are very well determined near halo masses of 
~ 1 X 10^'' h~^MQ and a redshift z ~ 0.2. The shght 
differences in these values for each experiment reflect the 
fact that the statistical errors have a different scale de- 
pendence for each experiment. LSST has the greatest sky 
coverage, so the LSST pivot mass is the largest. Con- 
versely, SNAP has much less sky coverage but with a 
galaxy number density of n = 100 arcmin^^ it is very 
deep and has the smallest mpiv. These mpiv and Zn 



^piv 



values are not unexpected. It is already well known that 
these masses are the largest contributors to the conver- 
gence power near 

The concentration constraints achieved in the course 
of this self calibration are stringent and potentially very 
useful. Marginalizing over the other parameters, LSST 
alone can provide a ~ 5% constraint on Cpiv at la, and 
can constrain the mass and redshift power-law indices to 
a{a) = 0.08 and cr(/3) = 0.27 respectively Both SNAP 
and DES constrain Cpiv at the ^ 10% level. 

Figure [5] shows concentration parameter constraints 
on the binned c(m, z) relation model. First, recall that 
our convention is to designate our concentration param- 
eters with a subscript of the value of the logarithm of 
the halo mass at the bin center. The parameter that 
describes halo concentrations at m = lO^^'^ H^^Mq is 
C13.5 and so on. Rather than displaying all parameter 
constraints. Figure [5] shows a reduced set of parameters 
that are at least mildly constrained by the experiments 
that we consider. The other concentration parameters 
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FIG. 7: Constraints on the effective halo concentration pa- 
rameters from tomographic weak lensing after marginalizing 
over the uncertainty in cosmological parameters. This plot 
shows Icr confidence contours on the parameters of the power- 
law concentration relation. We show the fractional uncer- 
tainty in Cpiv relative to the fiducial model. Recall that the 
pivot concentration is the concentration at the halo mass and 
redshift at which concentrations are best determined. In all 
panels, from outermost to innermost, the contours correspond 
to those achievable with DES, SNAP, and LSST using mul- 
tipoles up to ^max ~ 3000. In the upper right portion of the 
plot, we list the parameterization and the values of the pivot 
masses and redshifts for convenience. 



are constrained at uninteresting levels in all cases. 

The fact that halo concentrations are most well con- 
strained near halo masses near ~ 10^^ /i~^Mq can be 
seen directly in Figure [S] This is clearly the most well 
constrained parameter. In fact, this parameter can be 
constrained with an uncertainty cr(ci4.5)/ci4.5 = 0.48 rel- 
ative to the fiducial value by LSST alone. In the Ma 
et al. [i^ fiducial model with greater small-scale power, 
LSST constrains this parameter relatively more strin- 
gently with cr(ci4.5)/ci4.5 = 0.31. Generally, this param- 
eter is somewhat degenerate with the concentration value 
at the next lowest mass bin C13.5, and a combination of 
these two parameters is constrained at interesting levels. 
As with the power-law case, the redshift dependence is 
comparably poorly constrained. 

Though the constraints presented in Fig. [8] are mod- 
est, those in the power-law model of Fig. [7] are encourag- 
ing and recall that we have been rather conservative in 
our priors. Given that precise concentration constraints 
may inform the modeling of baryonic physics, it is inter- 
esting to estimate the best possible constraints on con- 
centrations that may be achieved with these methods. 
To minimize the interplay between dark energy and con- 
centration parameters, we may assume that dark energy 
parameters are known extremely well from other experi- 
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FIG. 8: Constraints on the effective halo concentration pa- 
rameters from tomographic weak lensing after marginalizing 
over the uncertainty in cosmological parameters. We show the 
usefully constrained combination of parameters in the C13.5- 
C14.5 plane and projections of these parameters with the red- 
shift variable /3. In all panels, from outermost to innermost, 
the contours correspond to the la contours achievable with 
DES, SNAP, and LSST with £m^^ = 3000. 
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FIG, 9: Same as Figure [71 but assuming prior constraints 
on cosmological parameters at levels that are expected from 
a combination of the Planck cosmic microwave background 
satellite and supernovae luminosity distance measures from 
SNAP, The additional prior matrix for this analysis is from 
Ref, [H, 



ments, for example through the supernovae measured by 
SNAP itself, and that we aim to measure effective halo 
concentrations via weak lensing. This is analogous to the 
dark energy constraints shown in Fig, [21 where it was as- 
sumed that concentrations were known perfectly, but in 
this case it is the cosmology that we assume to be well 
constrained. 

To estimate concentration constraints in this scenario 
realistically, we assume that we have prior knowledge of 
cosmological parameters that exceed the contemporary 
priors that we have thus far assumed. To be specific, 
we add priors that are characteristic of projected con- 
straints from the Planck cosmic microwave background 
mission and the supernova luminosity distance compo- 
nent of SNAP, In this case, the prior matrix does not 
have the simple diagonal form of our previous prior ma- 
trices. We have taken the prior Fisher matrices from the 
analysis of SNAP and Planck constraints in [s^ . 

We show constraints on effective halo concentrations 
in the power-law concentration model assuming much 
more restrictive priors on cosmological parameters in 
Figure [HI First, note the slight changes in the pivot 
masses and pivot redshifts. More importantly, reducing 
the dark energy degrees of freedom significantly tightens 
constraints on effective halo concentrations. For exam- 
ple, the la LSST constraint on Cpiv drops to about ~ 2%. 
Meanwhile, the LSST constraints on the mass dependent 
power-law index becomes cr(a) ~ 0.04. Restricting the 
space of cosmological parameters results in marginal im- 
provement on the redshift index, with a{(3) ~ 0.24. Both 



SNAP and DES also constrain concentrations stringently 
near the pivot mass and redshift, giving la constraints 
at the ^ 5% and ~ 7% levels respectively. 

The improvements with better cosmology priors are 
also noteworthy in the case of the general, binned c(m, z) 
relation. We show constraints on the subset of most well 
constrained parameters in Figure 1101 Consider the case 
of the LSST experiment, for which the constraints are 
the most stringent. In this case, the most significant 
improvement is realized for the redshift index /3, with 
a{(3) — 0.3, while the constraint on C13.5 decreases to 
just under 50%. The constraint on the concentration pa- 
rameter C14.5 from LSST drops to rather meaningful lev- 
els with cr(ci4.5)/ci4.5 ~ 0.3. Constraints in this regime 
are particularly sensitive to the true cosmology, which 
sets the signal-to-noise level of the forthcoming measure- 
ments. In the Ma et al. [i^l fiducial model with ag ~ 0.9, 
the parameters C14.5 and C13.5 can be independently con- 
strained to within 21% and 34% respectively with LSST. 



E. Constraint Degradation with Prior Knowledge 
of Halo Concentrations 

The experiments we have considered will report results 
beginning in the year 2011 in the earliest case of the DES 
or later. Of course, it is hard to predict the pace of the- 
oretical and observational studies of halo concentrations, 
but it seems reasonable to assume that over the next 
decade there will be some, perhaps significant, improve- 
ment and that this progress will be able to inform the 
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FIG. 10: Same as Figure |8l but assuming assuming prior 
constraints on cosmological parameters at levels expected 
from a combination of the Planck cosmic microwave back- 
ground satellite and supernovae luminosity distance measures 
from SNAP. The additional priors for this analysis are from 
Ref. i83|]. 



actual analysis of future weak lensing data. As such, it is 
interesting to explore the way in which prior constraints 
on the c(m, z) relation influence dark energy parameter 
degradation from tomographic weak lensing. 

In this section, we consider the influence of priors in 
the case of the binned c(m, z) relation. For simplicity, 
we assume that the concentrations in all mass bins, cn.s 
through C15.5, have a prior constraint of some fraction of 
their fiducial value. Let us designate this fraction dc- We 
illustrate the effects of prior knowledge of the c(to, z) rela- 
tion on parameter constraints by computing the parame- 
ter uncertainties using all information out to a maximum 
multipole £max = 3000 as a function of prior i5c. Though 
we take a single parameter to describe the prior in the 
interest of simplicity, it is clearly the prior on the pa- 
rameter C14.5 that is most important (see Sec. IIII D|) . We 
present results scaled relative to the uncertainties that 
would be derived under the assumption of perfect knowl- 
edge of effective halo concentrations, a^^'^''-^ . 

The dark energy parameter constraints with prior 
knowledge of effective concentrations is shown in Fig- 
ure [TT] As before, LSST is the most sensitive to concen- 
tration variations and can reap significant benefit from 
priors on the c(m, z) relation. In the LSST case, prior 
constraints on the concentration relation at redshift z = 
have little effect if they are greater than the ~ 30% level. 
This is sensible because LSST itself would constrain con- 
centrations in the relevant mass range at this level in the 
absence of prior knowledge. Naturally, prior constraints 
on c(m, z) that are better than ~ 30% lead to rapid de- 
creases in parameter uncertainties with decreasing 6c in 
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FIG. 11: The dependence of dark energy constraints on pri- 
ors on the c(m, z) relation. We plot the uncertainties in each 
parameter, scaled by the uncertainty under the assumption of 
perfect knowledge of the concentrations of halos a^^^''^ , as 
a function of the fractional prior on the amplitude of concen- 
trations Sc. As in Fig. [21 the rows from top to bottom show 
constraints on SIde, wq, and w^, while the bottom panels 
show the dark energy figure of merit, . The leftmost col- 
umn shows constraints with no prior on /3, the middle column 
shows constraints with a 20% prior on (3, and the rightmost 
column shows constraints with a 10% prior on /3. The line 
types are as in Fig. (2] namely, the solid line gives results for 
LSST, the rfas/ierf line corresponds to SNAP weak lensing, and 
the dotted line corresponds to DES. 



each case. 

The behavior of parameter constraints as a function of 
prior refiects the details of each of the experiments. An 
experiment with enormous sky coverage, such as LSST, 
can establish the amplitude of large-scale density fluctu- 
ations and the cosmic distance scale with low-multipole 
information and information from the high-redshift to- 
mographic bins rather effectively, so that such an exper- 
iment can take advantage of decreases in Sc even with 
no prior on /3. This statement is less true for SNAP 
and DES because these experiments have a relatively 
shorter lever arm in multipoles over which to derive con- 
straints. SNAP is limited by its comparably small sky 
coverage and DES by its comparable shallowness. As the 
prior constraints on /? improve, both DES and SNAP be- 
come increasingly capable of exploiting prior knowledge 
of the amplitude of the c(m, z) relation as can be seen 
for Sc ^ 0.25 in the two rightmost columns of Fig. [TT] 
In the case of LSST, the figure of merit increases with 
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decreasing 6c less rapidly than one might suppose based 
solely on the marginalized uncertainties in wq and 
because in this case the correlation between parameters 
varies strongly with Se- 
lla summary, the dark energy constraints from an ex- 
periment such as LSST can benefit greatly h'om prior 
knowledge of the effective concentrations of halos that 
are better than ^ 25%, while constraining the redshift 
evolution parameter /3, of concentrations to better than 
^ 30% leads to a rapid decrease in constraints in the ab- 
sence of prior knowledge of the amplitude of the c(to, z) 
relation. For DES and SNAP, it is necessary to constrain 
the redshift evolution of c(m, z) before either of these ex- 
periments can exploit prior knowledge of the amplitude 
of the c(m, z) relation effectively. In the cases we study, 
this requires prior knowledge of the P parameter to better 
than ^ 20%. Knowing the amplitude of the concentra- 
tion relation would reduce parameter degradation to less 
than ~ 20% in the case of LSST, while knowing both 
the amplitude and the redshift parameter (3 to less than 
~ 20% would reduce parameter degradation to less than 
~ 10% for all experiments. 



IV. SUMMARY AND DISCUSSION 

We have studied the ability of forthcoming weak lens- 
ing surveys to constrain both dark energy properties 
and the concentrations of dark matter halos simultane- 
ously. The primary motivations for this work are recent 
studies demonstrating that baryonic physics, neglected 
in prior forecasting of weak lensing constraints on dark 
energy properties, have a considerable influence on pre- 
dictions for convergence power spectra [13, HI]. More- 
over, the large enhancements in convergence spectra pre- 
dicted by these hydrodynamic simulations arose primar- 
ily from modifications in the internal structures of halos 
[28| . Clearly, gravitational lensing is sensitive to all mat- 
ter so here, as in the rest of the paper, we refer to halo 
concentrations, but by this we mean the effective concen- 
trations of composite systems comprised of dark matter 
and baryons in radial range 0.04 < f/_R200m ^ 1- We 
have shown that neglecting these baryonic effects can 
lead to large biases in inferred dark energy parameters 
unless information from multipoles greater than a few 
hundred is disregarded (see Fig. [3]). On the other hand, 
weak lensing surveys have the ability to self calibrate for 
the effects baryons as an alternative to eliminating small- 
scale information. 

We modeled the unknown relationship between effec- 
tive concentration and halo mass in two different ways. 
In the more restrictive set of models, we considered con- 
centration as a power-law function of halo mass. As we 
have demonstrated, convergence spectra are most sensi- 
tive to halo concentrations in a relatively narrow range 
of halo masses near a few xlO^^ h~^MQ, so this class 
of models may not be overly restrictive. Nevertheless, in 
currently-favored cosmological models, the mass-to-light 



ratios of halo systems must be a non-monotonic func- 
tion of halo mass, so it is desirable to study relations 
with more freedom. Along these lines, we studied binned 
relationships between concentration and halo mass spec- 
ified by average concentration values at a fixed number 
of halo masses, evenly spaced in the logarithm of the 
halo mass between 11 < log(A'// h~^MQ) < 16. In these 
models, we used spline interpolation to obtain halo con- 
centrations between the values of halo mass at which the 
concentrations were specified. Throughout this paper, 
we have presented results for a concentration-mass rela- 
tionship specified at five bins in this mass interval, but 
have verified that our results are relatively insensitive 
to the number of bins up to about twelve bins. Coupled 
with our prejudice that these effective concentrations will 
likely be very smooth functions of halo mass, this class of 
model is likely to have sufficient fiexibility to produce the 
types of mass-concentration relations that are realized in 
nature. In both cases, we allowed for a redshift depen- 
dence to the modification in the halo mass-concentration 
relation by taking c(m, z) = c{m, z — 0)/(l -I- z)^ . 

Our study implies that self-calibration is a quite 
promising alternative to simply eliminating multipoles 
beyond which significant parameter biases may be likely. 
Comparing constraints obtained by removing multipoles 
beyond a few hundred with constraints that include all 
multipoles up to £ = 3000 and assume perfect knowledge 
of halo concentrations, the degradation in parameter un- 
certainties is a factor of ~ 1.5 — 3. The details of the 
degradation depend upon the parameter of interest, the 
experiment, and the fiducial cosmological model. As it 
is the most sensitive of the forthcoming instruments on 
large scales, the requirements for LSST to be unbiased are 
the most restrictive and indicate that multipoles beyond 
^ 300 should not be considered in parameter extraction. 
The corresponding degradation in parameter constraints 
for LSST is a factor of '--^ 3 for all parameters. On the con- 
trary, our study suggests that self calibration is a more 
fruitful alternative. Even in our most permissive concen- 
tration models, dark energy parameter constraints from 
this self-calibration exercise are never degraded by more 
than 36% and more restrictive models (if they can be 
justified) do significantly better. In particular, if the 
concentration-mass relation can be constrained to better 
than 20% at a specific redshift, the constraints derived 
from an LSST-like experiment will be degraded by only 
^ 20%. If, in addition, the power-law index describing 
the redshift evolution of effective halo concentrations is 
also known to within 20%, then dark energy parameter 
constraints will be degraded by no more than ~ 10% for 
any of LSST, DES, or SNAP individually. 

Self calibration as an alternative to removing the in- 
formation contained in high-£ multipoles has another ad- 
vantage; it provides direct constraints on halo concentra- 
tions that may complement, if not be competitive with, 
other methods. In the case of the power-law class of 
models for effective halo concentrations as a function of 
mass and redshift, LSST alone can constrain halo con- 
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centrations near the pivot mass, rripiv ~ 10^^ h~^MQ, 
and redshift, Zpiv ~ 0.2, at the ~ 5% level with only 
modest prior constraints on cosmology that are typical 
of our contemporary knowledge of cosmological parame- 
ters. Both DES and SNAP constrain halo concentrations 
at the pivot mass and redshift to about ~ 10%. General- 
izing to the case of the arbitrary, binned relation for halo 
concentrations, LSST alone constrains halo concentra- 
tions near to 3 x 10^"* /I'^M© to within 48%. Adding 
prior constraints on cosmological parameters from the 
Planck cosmic microwave background satellite and su- 
pernova luminosity distances as measured from SNAP 
reduce this to a 30% constraint. Though modest, this 
constraint would not be subject to the same systematics 
and selection issues as other techniques. 

Our work has several shortcomings. Of course, our 
framework for all of these calculations has been the halo 
model, and this phenomcnological model has known inad- 
equacies. More importantly, we have modeled the effect 
of baryons as entirely restricted to the internal struc- 
tures of halos. Thou gh t his is a good approximation 
in current simulations [28|, it is unreasonable to expect 
that this is strictly true and no study has yet addressed 
potentially-important, though sub-dominant, additional 
modification to matter power spectra. Changes to the 
detailed profile shapes, and the mass function of halos 
(already measured in Ref. the clustering of halos 

are all likely to be important at some level. To address 
this issue is an ambitious goal which lies well beyond 
the scope of our study, yet it is a goal that must be 
achieved in order to take full advantage of these forth- 
coming experiments. The degree to which baryons will 
modify the matter power spectra will likely remain un- 
certain for some time, yet we must at least enumerate the 
different ways in which the power spectra are modified. 
As an example, the study of Ref. [2^ does show percent- 
level residuals between their simulation results and the 
results of a model based entirely on modified halo con- 
centrations. Although the results of Ref. 28] are noisy 
due to their small simulation volume, we can use their 
data to obtain a first guess at the residual bias that may 
remain after allowing for variation in halo concentration 
parameters. Our estimates suggest that in the case of 
LSST the bias is a few percent of the la parameter un- 
certainty at £max — 10'^ and is of order \S\/a ~ 0.1 at 
^max = 3000 for both wq and Wa- The residual biases in 
the case of DES and SNAP are smaller yet. Furthermore, 
we have described forthcoming experiments in a simple 
manner and so detailed comparisons between the pro- 
jected performances of instruments based on this study 
should not be regarded as definitive. Certainly, system- 
atic issues regarding shape measurements will be better 
controlled in space-borne instruments like SNAP than 
from the ground. Moreover, SNAP with potentially much 
greater depth than either DES or LSST will be sensitive 
to relatively higher-redshift, larger-scale fluctuations at 
a fixed angular scale. On such scales, the fiuctuation 
spectrum will be less subject to nonlinear effects. Lastly, 



we have ignored all additional systematics whether theo- 
retical or observational. We anticipate that in the years 
leading up to the surveys we have studied, a considerable 
amount of effort will be put into making theoretical pre- 
dictions more accurate. We expect — or at least hope — 
that more sophisticated methods for predicting matter 
power spectra will supplant current techniques and find 
it most likely that the parameterizations we adopted will 
not be those that are implemented in the actual data 
analyses. Nevertheless, our study gives a preliminary in- 
dication that calibration of uncertain baryonic effects will 
be possible at minimal cost and that this calibration may 
provide interesting information about halo structure that 
may be used to diagnose models of the evolution of the 
baryons themselves. 

In the context of our models, calibrating effective halo 
concentrations in parallel with deriving constraints on 
dark energy parameters is extremely successful and the 
preceding paragraph is not intended to undercut our pri- 
mary results and conclusions. We are optimistic that 
this success will carry over to more complex models that 
may include systematics and additional modifications to 
power spectra due to baryonic processes. Part of our op- 
timism stems directly from the success with which halo 
structure can be calibrated, as we have demonstrated 
here. In addition, we are optimistic because we have 
used only a fraction of the information contained within 
these surveys. In particular, we restricted our analysis 
to multipoles £ < 3000 in order to remain in the regime 
where Gaussian statistics and other weak lensing approx- 
imations are valid. Accounting for non-Gaussianity and 
dropping additional assumptions, such as the Born ap- 
proximation and the use of reduced shear (66j . would 
allow higher multipoles to be included in the analysis. 
This would greatly increase the constraining power of 
these week lensing surveys. More specifically, SNAP is 
the deepest survey we have considered by a large margin. 
As such, SNAP suffers the most from setting ^max — 3000 
and going beyond this limit would enable SNAP to out- 
perform the projections we present here (while system- 
atics will also be better controlled from space). An im- 
portant implication of this is that upon extending the 
analysis to small scales, SNAP would provide the best 
determination of baryonic effects, eliminating degenera- 
cies that may limit dark energy constraints. Of course, 
gains from small-scale observations would also come with 
the challenge of predicting power spectra out to higher 
multipoles. Furthermore, we have used only weak lensing 
shear information and have not utilized the galaxy power 
spectra and galaxy-shear cross spectra that would also 
be contained within these surveys and would make these 
constraints more robust [6^ [sj] • Our study shows the 
great promise of future surveys to constrain dark energy 
along with numerous other phenomena, but is is only the 
tip of the iceberg and we must meet several theoretical 
challenges in order to exploit the enormous amount of 
information that will be available in forthcoming photo- 
metric surveys. 
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